home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / os2 / adaptor.zip / ADAPT.ZIP / adaptor / examples / simple / lprime.f < prev    next >
Text File  |  1993-06-28  |  794b  |  35 lines
  1.       program prime
  2.       integer n, s, k
  3.       parameter (n=10000)
  4.       logical*1 a(n)
  5. !hpf$ distribute a(block)
  6.       call cm_timer_clear (1)
  7.       call cm_timer_start (1)
  8. !hpf$ independent, local_access
  9.       do i = 2, n
  10.         a(i) = .true.
  11.       end do
  12.       a(1) = .false.
  13.       k = 2
  14.       do while (k*k .le. n)
  15. !hpf$    independent, local_access
  16.          do i = k*k,n,k
  17. c          all multiples of k are no primes
  18.            a(i) = .false.
  19.          end do
  20.          k = k + 1
  21.          do while (.not. a(k))
  22.            k = k + 1
  23.          end do
  24.       end do
  25.       s = 0
  26. !hpf$ independent, local_access
  27.       do i = 1, n 
  28.         reduce (count, s, a(i))
  29.       end do
  30.       call cm_timer_stop (1)
  31.       print *, 'There are ',s,' primes until ', n
  32.       call cm_timer_print (1)
  33.       end
  34.  
  35.